ONE UPPER ESTIMATE 
ON THE NUMBER OF LIMIT CYCLES 
OF EVEN DEGREE LIENARD EQUATIONS 
IN THE FOCUS CASE 

GRISHA KOLUTSKY 

Abstract. We give an explicit upper bound for a number of limit 
cycles of the Lienard equation x = y — F(x), y = —x of even degree 
in the case its unique singular point (0,0) is a focus. 

M. Caubergh and F. Dumortier get an explicit linear upper 
estimate for the number of large amplitude limit cycles of such 
equations [CD . We estimate the number of mid amplitude limit 
cycles of Lienard equations using the Growth-and-Zeros theorem 
proved by Ilyashenko and Yakovenko [IYa . 

Our estimate depends on four parameters: n, C, a±, R. Let 

n-1 

F(x) = x n + diX 1 , where n = 21 is the even degree of the monic 

polynomial F without a constant term, Vi \a,i\ < C, so C is the 
size of a compact subset in the space of parameters, R is the size of 
the neighborhood of the origin, such that there are no bigger than 
I limit cycles located outside of this neighborhood, |oi| stands the 
distance from the equation linearization to the center case in the 
space of parameters and 2 — \a\ | stands the distance from the equa- 
tion linearization to the node case in the space of parameters. 
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I. Hilbert-Smale problem 

In 1977 A. Lins Neto, W. de Melo and C. C. Pugh jLMPj examined 
small perturbations of a linear center for a special class of polynomial 
vector fields on the plane. This class is called Lienard equations: 

(x = y-F(x), 
{ y = -x, 

where F is a polynomial of odd degree. Actually, Lienard in 1928 in- 
troduced it for a modeling of the non-linear damping in electric circuits 
[L] . It was a generalization of the famous Van der Pol equation [V] . 

Authors of |LMPj proved the finiteness of limit cycles for a Lienard 
equation of odd degree n. Let us remind that the Finiteness problem 
(also known as the "Dulac problem") was solved in full generality only 
in 1991 by Ilyashenko [IT] and in 1992 by Ecalle [E] independently. 

Also A. Lins Neto, W. de Melo and C. C. Pugh jLMPj conjectured 
that the number of limit cycles of ([l| is not bigger than 

In 1998 S. Smale [S] suggested to consider a restriction of the second 
part of the Hilbert's 16th problem to Lienard equations of odd degree. 
He conjectured that there exists an integer n and real C such that the 
number of limit cycles of ([!]) is not bigger than Cn q . 

In 1999 Yu. Ilyashenko and A. Panov [LP] got an explicit upper bound 
for the number of limit cycles of Lienard equations through the (odd) 
degree of the monic polynomial F and magnitudes of its coefficients. 
Their result reclined on the theorem of Ilyashenko and Yakovenko that 
binds the number of zeros and the growth of a holomorphic function 

ra. 

In 2007 F. Dumortier, D. Panazzolo and R. Roussarie jDPRj con- 
structed a counterexample to the conjecture of A. Lins Neto, W. de Melo 
and C. C. Pugh. Namely, they presented an example of a Lienard equa- 
tion of odd degree n with at least limit cycles. 

In 2008 Yu. Ilyashenko [13] suggested to prove a result analogous to 
the one of Ilyashenko and Panov for Lienard equations of even degree. 

In 2008 M. Caubergh and F. Dumortier in jCD] proved the following 
theorem for Lienard equations of even degree. 

Theorem 1. Let K be a compact set of polynomials of degree exactly 
n = 21, then there exists R > such that any system having an expres- 
sion (IT]) with F e K has at most I limit cycles having an intersection 
with W\B R . 

Here and bellow Br denotes the ball around the origin with the 
radius R. 

2. Notations and the Ilyashenko strategy 

From now on we will consider a system 0, where F is a monic 
polynomial of even degree n = 21 without a constant term. 
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Remark 1. The assumption F(0) = does not reduce the generality; 
it may be fulfilled by a shift y i— > y + a. The assumption that F is monic 
may be fulfilled by rescaling in x, y and reversing the time if necessary. 

Let v be an analytic vector field in the real plane, that may be 
extended to C 2 . For any set D in a metric space denote by U £ (D) the 
e- neighborhood of D. The metrics in C and C 2 are given by: 

p(z,w) = \z — it? I, z,w E C; 

p(z,w) = maxd^i — wi\, \z2 — W2I), z,w E C 2 . 

Denote by \D\ the length of the segment D. For any larger segment 
D' D D, let p(D, dD') be the Hausdorff distance between D and dD'. 

We want to apply the next theorem proved by Ilyashenko and Panov 
|IPj . In fact, it is the easy corollary from the Growth-and-Zeros the- 
orem for holomorphic functions proved by Ilyashenko and Yakovenko 

Consider the system 

x = v(x), xER 2 . (2) 

Theorem 2. Let T be a cross-section of the vector field v, D C T 
a segment. Let P be the Poincare map of ([2| defined on D, and 
D C D' = P(D). Suppose that P may be analytically extended to 
U = U £ (D) c C, e < 1, and P(U) C U\D') C C. Then the number 
#LC(D) of limit cycles that cross D admits an upper estimate: 

The same is true for P replaced by P -1 . 

Actually, the Ilyashenko strategy is the application of the previous 
theorem. It requires a purely qualitative investigation of a vector field, 
i.e. a construction of such D for every nest of limit cycles. This strategy 
was applied before in papers [12] and [TP] . 

We take K from the Theorem [T] to be the space of monic polynomials 
of degree exactly n with coefficients, which absolute values are bounded 
by some positive constant C > 4, i.e. 

n-l 

F(x) = x n + y^diX 1 , Vz : \ai\ < C. 

i=l 

If < |ai| < 2 then the unique singular point (0, 0) of the system (JTJ) 
is a focus. In our work we will consider only this case. 

3. Bendfxson trap from within 

In this Section we construct an interval D, which lies inside Br and 
intersects transversally all limit cycles in Br. Also we find an upper 
estimate for the Bernstein index, b = log (n an/) . To do that we need 
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to estimate p(D, dD') from bellow, where D' = P(D) C D and P is the 
Poincare map defined on D (see the Figure [T]) . 




Figure 1. The inverse Poincare map of the Lienard 
equation ([I]) inside the ball B R . 

Let (p, r be polar coordinates on M 2 , ip, r be derivatives with respect 
to®. 

First of all we need to determine the size of the domain, there the 
Poincare map is defined. 

Lemma 1. Put a = l ai l(jH ai D e iaii-2 _ j n tfo e f ocus case (o < \a x \ < 2) 
the Poincare map for the system ([I]) is well defined in B a . 

Proof. Let us calculate f. 

xx + yy r cos tp{r snap — F(r cos tp)) — r 2 sin ip cos (p 



cos ipF(r cos ip>) = — r cos 2 ip Oj(r cos ip) 



i-l 



-r cos 2 ip (at + 0(r, <p)) , (4) 



where 0(r,<p) = a i( r cos f) 

i=2 

Let us calculate (p. 



i-l 



xy — yx —r 2 cos ip — r 2 sin 2 <p + r sin <pF(r cos </?)) 



-1 + 



sin (pF(r cosy?) 



-l +si n2 H | + ^). (5) 
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The absolute value of the function 0(r, ip) admits the following upper 
estimate in Bi: 

2 

" i _ r «-i n r 

\0(r, <p)\ < V Cr 1 - 1 = Cr < < 2Cr. (6) 

u 1 - r l-r ~ v/ 

1=2 

Therefore, in B 2 -\ ai \ cBi: < l ai j~ 2 . Indeed, 

AC 2 

^<-l + |sin2 P | (M + M^T)<_ 1 + M +Cr < 

loil — 2 2 — lail lad — 2 
< J — ^ h C- 



4C 



Also, in B 2 -\ ai \ : |r| < 2r. Indeed, 

4C 



|f| < |rcosVl Oil + |0(r,y)|) < (K| + 2Cr)r < 

f, ^2 - |ai|\ 2 + lad 

< (\ ai \+2C—^Jr<—^-r<2r. 

Hence, any trajectory starting from any point from B a rotates around 
the origin on the angle not less than 27r before leaving -B g-iajj . 

Indeed, (p < l ai j~ 2 implies that during the time, A t = 2ir 2 _ 4 ^ the 
variation of the angle, A ip > 2tt and the variation of the radius, 
A r < e 2_ l a ii during the same time A t, because |r| < 2r. 

Finally, ae^\ = JM|dM < ?=N. □ 

Let us denote by Y" the maximal ^/-coordinate of the point of inter- 
section between the most external limit cycle which lies inside Br (if 
it exists, of course) and y-axis. 

Lemma 2. If ai is negative, then r > in B a . Let D = [cr,Y] C Oy. 
Thend = p(D,dD') > ^a. 

Proof. If r < cr, then r < \ and by Q: |0(r,p)| < 2Cr < ^. 
Therefore by Q, 



r > rcos <p I — a-y —J = — — rcos > 0. 

This proves the first part of the Lemma. 

Consider the orbit 7 of the system Q that passes through the point 
(0, cr). Then the Hausdorff distance, d can be estimated as follows: 



d > 



2- 



f^dcp 



2tt 



> / — — a cos ipdip = — - — cr. 







This inequality completes the proof of the Lemma. □ 

Remark 2. For positive a\ we can get the same results just by reversing 
of the time. 
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Now we can estimate b from above: 

R + 2{R + 2) R + 2 

b<\og — :— <log^ — < -: — j— . (7) 
a TT\ai\a \a>i\cr 

4. Complex domain of the inverse Poincare map 



The Theorem [2] uses the width e of the complex domain U e (D) to 
which the (inverse) Poincare map may be extended. We will apply the 
following theorem to estimate this e from bellow. 

Theorem 3. Let P : D — > D' be the Poincare map of (J2J). For any 

x E D denote by <p x ,p{x) the arc of the phase curve of ([2| starting at x 
and ending at P(x). 
Let 

Q(D) = [J cp x ,p( x ), 



and 



1 < a = max \v\, L = 2u. (8) 

~ c/2(n) 



Let t(x) be the time length of the arc tp x ,p{x), and 

^~max = maxtfx), T = T max + 1. 
xeD 

Let 

6<e~ LT , \ = V6, e = 5 2 . (9) 

Suppose that (^1,^2) are coordinates in C 2 , c r = {z\ = 0}, 
v = (v 1 ,v 2 ). 

Let U s = U 5 (0) x U X (D') C C 2 . Suppose that 



V2 



<H in U s . (10) 



Then the Poincare map P : D — > D' of ^ may be analytically ex- 
tended to U £ (D) C c r ; and P{U £ {D)) C U\D). 

The same is true for P replaced by P~ x . In this case P~ X {D) = D' , 

tt(D) = |J (p x ,P(x). 
xeD' 

For the proof see [IP]. □ 

Bellow we will produce some preliminary calculations, which would 
allow us to apply the Theorem [3] later. 

Definition 1. A C-monic polynomial is a real polynomial in one vari- 
able with the highest coefficient one and other coefficients no greater 
than C in absolute value, with zero constant term. 
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Proposition 1 (Properties of C-monic polynomials). Let F be a 

C-monic polynomial of degree n, C > 2. Then 

max \F(x)\ < 2X n for X > C + 1, (11) 

xe[o,x] 

max \F'(x)\ < Cn 2 X n ~ l for X > 1, (12) 

xE[0,X] 

\F(z)\<2C\z\ for z G C, \z\ < -. (13) 

For the proof see [IP] . □ 



Lemma 3. Let v be the vector field given by the system (fTl). Then jjt, 
and L from the Theorem [3] admits the following estimates: 

(i<3{R + 2) n L<6(R + 2) n . (14) 

Proof. By definition, U 2 {Vl) C B R+2 . So 

\v\ < \x\ + \y\ < \x\ + \y\ + \F(x)\ < 2(R + 2) + 2{R + 2) n , 



■ 



where the last inequality provided by (11). Hence, 

fi<2(R + 2 + (R + 2) n ) < 3(R + 2) n , L = 2fi < 6(R + 2) 
that proves the Lemma. □ 
Let G = B R \B a . Then O = |J (p X)P ( x ) C G. 

Lemma 4. Let ^ y be the arc fy,p-i{ y ) of the phase curve of where 
y G D. Then t(y), the time length ofj y , admits an estimate 

, x 25C 2 n 2 R n , . 

T max = maxt(y) < . (15) 

Proof. The arcs j y , y G D belongs to G. We will split G into two 
domains: |x| < a and \x\ > a for a small to be chosen later. The 
second domain contains two parts of 7 y : one with x < —a, the other 
with x > a. The time length of any of them is no greater than — . 
In the next Proposition we will choose a so small that the curvilinear 
strip 

S a = {(x,y) EG: \y - F(x)\ < a} 
is crossed by the orbits of ([TJ in the time no greater than 1. 

Proposition 2. Let 

U 3C" a 2Cn 2 R n ~ 1 QC 2 n 2 R n - 1 ' 1 ' 

Then the time length of any arc of the orbit of Q located in S a is no 
greater than 1. 
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Proof. By the symmetry arguments it is sufficient to prove that in 

s+ = s a n {x > 0}: 

j t (y - F(x)) < -2a. 

Let us first prove that in S a we have: \x\ > uj. Namely, let |x| < uj, 
\y — F(x)\ < a. Then (x,y) E D a . Indeed, 

\x\ + \y\ < uj + a + max 

[0,ui] 



By (jTej) , a < u < \. By (g, \F(x)\ < 2Cuj. Hence, for x e [0,^], 
|ar| + \y\ < (2C + 2)w < 3Cu; = a. 

By g, < Cn 2 ^" 1 in G. Therefore, for x such that (a;, y) e S+ 

we have: x > uj, and 

d 



dt 

because 



(y - F(x)) = -x- F'(x)(y - F(x)) < -uj + aCn 2 ^"" 1 < -2a. 

UJ UJ 



a 



< 



2Cn 2 R n ~ 1 2Cn 2 R n ~ 1 + 2 



□ 



Let us finish the proof of the Lemma |4j 

The arc 7 y spends in S a no longer time than 2 (two crossings, each 
one no longer in time than 1, by the previous Proposition); in G \ S a 
no longer time than — (two crossings, one to the left, another to the 
right with \x\ > a). Hence, 

AR ■ QC 2 n 2 R n ~ l 25C 2 n 2 R n 

T < 9 A < 

± max _i ^ i ~ 



a 



a 



This calculation completes the proof of the Lemma |4| □ 

Remark 3. The same inequality holds for T max replaced by T max + 1. 

Let us check the last assumption of the Theorem [3j 

Lemma 5. Take 

300C 2 n 2 R n (RAr2y 



cxp 



a 



, S = y/e, \ = VS. (17) 



Let, as in the Theorem^ Ii & = U\0) x U (D') C C 2 . Then in n^: 

v 2 

— < ft. 

Vl 

Proof. By (|13j) and by definition of n,5, 

\vi(z)\ > \\z 1 \-F(z)\ > (a- A) -2CS > 6, 
where the last inequality is trivial. On the other hand, v% : 



n«5, \v 2 \ < 5. Hence, 



< 1 < jJL. 



-x. In 

□ 
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Lemma 6. The inverse Poincare map of the Lienard equation ([I]) may 

_ 300C 2 n 2 R n (R+2) n 

be extended to the domain U £ (D) C C, where e = e 
Moreover, p-\U e {D)) C U l {D'). 

Proof. This Lemma following from the Theorem [3j Lemmas [3] and [5] 
verifies assumptions [8] and [10] respectively. We only should check the 
assumption 9j By the Remark |ij T < 25C ' 2 ™ 2Rn . Hence, 



cxp 



150C7 2 n 2 i? n ( J R + 2) n 



/ 25L7 2 n 2 iT\ . _ L , ; 



exp -6(i2 + 2) n < e 



(7 



that proves the Lemma. □ 



5. Final estimate 

Theorem 4. The number L(n,C,ai, R) of limit cycles of (fl| in the 

case when n is even, C > 4 and < |ai| < 2, admits the following 
upper bound: 

T . „ t~*\ ( /38400C 4 n 2 J R n+1 ( J R + 2) n+1 ^ 
L{n, C, ai, i?) < exp exp . — — . — — e 2 -i«i 

V V l°il ( 2 - Kir 



Proof. Now we can apply the Theorem |2| By definition, \D\ and \D'\ 
are less than R. The Lemma [6] provides us with the lower bound on e. 
So estimates ^ and ^ imply: 



Ti „ ( r, {300C 2 n 2 R n (R + 2) n \\ R + 2 
L(n,C,a!,R) < exp I 2Rexp I — ^ '— \ J 

(2R(R + 2) r300C 2 n 2 R n (R + 2) n \\ 
< exp — p — : exp — < 

V l a il°" V ° ) ) 

( f600C 2 n 2 R n+1 (R + 2) n+1 
< exp ' 



< 



I exp I - 



ladcr 2 



'38400C 4 n 2 J R n+1 (i? + 2) n+1 / 16tt 
exp l exp l : — r— 7- ; — rr^ exp 



|«i | 3 (2 — |) 2 ^V2-h| 
This calculation completes the proof of the Theorem. □ 
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